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PREFACE 

This  report  is  submitted  in  compliance  with  the  requirements  of  the  U.S. 
Air  Force  Special  Weapons  Center  Contract  Number  F2960i-7^C-0055  "Lightweight 
Alternator"  by  the  AiResearch  Manufacturing  Company  of  California,  a division 
of  The  Garrett  Corporation,  Torrance,  California. 

The  report  discusses  the  test  and  computer  simulation  activities  conducted 
on  a mod i f ied-des i gn  , small-rating  aircraft  alternator.  These  activities  were 
intended  to  provide  additional  verifications  of  the  alternator  analytical 
approach  utilized  in  designing  two  specialized  large-rating  lightweight  alter- 
nators. The  latter  designs  have  been  previously  detailed  in  Air  Force  technical 
reports  AFWL-TR- 75-66 , "Subsystem  Design  Analysis  Report  for  Lightweight 
Alternator  (AC  Load  Case),"  and  AFWL-TR-75-66 , Addendum  1,  "Subsystem  Design 

M ' 

Analysis  Report  for  Lightweight  Alternator  (DC  Load  Case)". 

Detailed  technical  direction  was  provided  by  the  Air  Force  Special  Weapons 
Center,  Kirtland  Air  Force  Base,  Albuquerque,  New  Mexico.  Major  Frank  Zimmer- 
mann  and  Mr.  James  O'Loughlin  provided  continuing  direction,  guidance,  and 
support.  The  program  manager  at  AiResearch  was  Mr.  T.E.  Brown.  Principal 
investigators  and  contributors  on  technical  activities  at  AiResearch  were: 

Dr.  C.H.  Lee,  Mr.  D.  Berker,  Mr.  G.  Tatro,  and  Mr.  P.  Walla. 
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SECTION  I 

OBJECTIVES  AND  APPROACH 

This  addendum  report  documents  the  final  activities  on  Contract 
F29601 -7^-C-0055  consisting  of  alternator  test  and  computer  modeling  activities 
accomplished  in  accordance  with:  (l)  the  Air  Force  Statement  of  Work  (Purchase 

Request  FY3592-75-10147,  Amendment  No.  2 to  Statement  of  Work  for  Purchase 
Request  FY3592-74-10084) , and  (2)  the  technical  approach  and  program  specified 
in  AiResearch  proposal  75-11176,  dated  January  28,  1975. 

The  primary  objective  of  this  work  was  to  demonstrate  and  verify  the 
alternator  analytical  methods  and  techniques  utilized  to  attain  the  full-scale, 
large-rating,  lightweight  alternator  designs  previously  conducted  and  described 
in  the  preceding  reports  on  Contract  F29601-74-C-0055. ^ 

A secondary  objective  of  this  work  was  to  demonstrate  the  ac  resonant 
loading  concept.  Meeting  these  two  objectives  involved  correlation  of  results 
from  an  experimental  (test)  activity  and  a companion  computer  modeling  and  I 

analytical  activity.  Achieving  thorough  quantitative  definition  of  the  alter-  j 

nator  that  was  tested  and  modeled,  and  realizing  complete  and  accurate 
dynamic  modeling  methods  that  simulated  actual  tests,  were  of  oaramount 
importance  in  meeting  the  primary  work  objective. 

The  technical  approach  used  to  accomplish  the  above  objectives  consisted 
of  the  following  major  tasks: 

(a)  Analyze  an  existing  small  rating  test  alternator;  determine  needed 

modifications  to  simulate  the  previously  designed  full-scale  ; 

alternator  hardware  configuration  and  operating  characteristics,  i 

; 

and,  as  the  final  step,  quant i tat ivel y define  the  modified 

alternator.  ’ 

(b)  Modify  the  small  rating  alternator  per  the  final  step  in  (a)  above  I 

and  construct  a special  loading  and  control  circuit.  I 


Reference  1;  AFWL-TR-75-66,  Subsystem  Des ign  Analysis  Report  for  Lightweight 

A1 ternator  (ac  Load  Case) 

Reference  2:  AP/JL-TR-75-66,  Addendum  1,  Subsystem  Design  Analysis  Report  for 

Lightweight  Alternator  (dc  Load  Case) 
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(c)  Test  the  modified  alternator  in  combination  with  the  special 
loading  circuit  under  specified  transient  load  conditions  to 
experimentally  determine  the  alternator  characteristics. 

(d)  Construct  a simulation  model  of  the  modified  alternator  and  the 
specialized  load  circuit  using  the  Air  Force  SCEPTRE  transient 
analysis  computer  program.  Quantitative  basis  for  the  alternator 
model  was  the  design  description  obtained  as  the  final  step  in 

(a)  above. 

(e)  Using  the  model  developed  in  (d)  above,  analytically  determine 
the  alternator  performance  under  the  specified  transient  load 
conditions  and  compare  these  results  with  the  actual  alternator 
test  results  obtained  from  (c)  above. 

The  small-scale  alternator  used  was  an  existing  and  available  6-kva, 
400-Hz,  12,000-rpm  (4-pole),  120/208-v,  self-excited  laminated  salient  pole 
rotor  aircraft  type  of  machine.  Basic  modifications  made  to  the  machine 
cons  i s ted  of : 

(a)  Constructing  a replacement  solid  nonsalient  pole  rotor  for  the 
machine,  closely  simulating  the  rotor  configuration  selected  for 
the  full-scale  machinery  previously  designed  under  Contract 
F29601-74-C-0055. 

(b)  Providing  direct  excitation  of  the  machine  by  removing  and  replac- 
ing its  self-excitation  provisions. 

(c)  Adding  to  the  machine  instrumentation  needed  to  monitor  key  test 
informat  ion. 

The  modified  machine  design  description  was  established  using  the 
AiResearch  digital  computer  design  program  ROUND,  which  also  was  used  to 
develop  the  previous  full-scale  alternator  designs.  The  replacement  rotor 
designed  and  constructed  was  a solid-core,  nonsalient  pole,  distributed-winding 
configuration  without  rotor  amortisseur  provisions  (i.e.,  the  same  basic 
configuration  selected  for  the  full-scale  designs).  The  design  constants 
developed  using  ROUND  were,  as  noted  above,  utilized  in  the  later  computer 
modeling  and  transient  behav  lor  simulation  activities. 

The  modified  alternator  first  was  separately  tested  for  basic  character- 
istics and  performance  using  procedures  similar  to  those  defined  in  the  IEEE 
Mo.  115  Test  Procedures  for  Synchronous  Machines.  This  testing  induced 
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determination  of  saturation  character  is t ics , design  point  verification,  and 
definition  of  key  machine  reactances  and  constants  for  comparison  with  those 
developed  for  the  design  by  the  ROUND  analytical  program. 

A specialized  loading  circuit  and  controls  were  designed,  constructed, 
checkout  tested,  and  then  utilized  in  combination  tests  with  the  modified 
alternator  to  obtain  actual  transient  performance  of  the  machine.  An  original 
objective  for  this  loading  circuit  was  to  demonstrate  the  ac  resonant  loading 
concept  at  precise  resonance  conditions.  This  specific  objective  was  not 
met;  however,  the  machine  was  tested  at  two  off- resonance  conditions. 

Two  different  transient  analysis  computer  programs  (i.e.,  SCEPTRE  and 
an  AiResearch  program)  were  utilized  to  perform  transient  simulations  of  the 
alternator  and  specialized  load.  These  two  programs  differ  in  solution 
methods  but  both  handle  the  same  basic  mathematical  model.  In  this  model, 
the  alternator  is  represented  by  a six-circuit  equivalent,  consisting  of  three 
armature,  one  field,  one  direct  camper,  and  one  quadrature  damper  circuit. 

The  latter  two  circuits  in  this  case  represent  the  damping  action  of  the 
solid  rotor  core  (i.e.,  no  discrete  amortisseur  or  damper  windings  are  included 
in  this  rotor). 

Treatment  of  the  SCR  switching  network  is  different  in  the  SCEPTRE  and 
AiResearch  analysis  programs.  In  the  SCEPTRE  program,  the  SCR's  are 
approximated  by  resistors  in  series  with  the  load  capacitor.  The  resistor 
value  is  appropriately  varied  between  two  set  extremes  so  that  it  has  a 
value  near  zero  ohms  when  the  capacitor  is  charging,  and  a high  value 
(e.g.,  10,000  ohms)  when  the  capacitor  is  maintaining  or  cischarging.  In 
the  AiResearch  program  the  SCR  is  directly  treated  as  a switch  that  may  be 
set  open  or  closed. 

Modeling  treatment  of  the  effect  of  the  magnetic  saturation  on  the  model's 
inductive  values  is  also  different  in  these  two  programs.  Inherently,  the 
SCEPTRE  program  requires  that  the  inductance  m.atrix  be  symmetrical.  This  is 
not  a prerequisite  in  the  AiResearch  program.  The  numerical  method  of 
solution  used  in  the  AiResearch  program  is  a miodified  or  improved  Euler's 
technique.  The  SCEPTRE  program  can  utilize  any  one  of  three  explicit  integrating 
techniques:  Runga-Kutta,  trapezoidal,  or  XPO  technique.  The  latter  has 

been  utilized  in  conducting  this  work. 
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Comparison  of  results  using  these  two  programs  shows  no  essential 
differences,  that  is,  both  programs  yield  identical  solutions.  The  Euler 
technique  utilized  in  the  AiResearch  program  results  in  a negligible  sacrifice 
in  precision  (i.e.,  compared  with  the  SCEPTRE  program  XPO  technique),  which  is 
more  than  offset  by  reduced  computational  and  compilation  time. 

Overall,  results  of  the  alternator  tests  and  computer  simulations  correlate 
very  well.  Load  cases  were  established  for  both  the  test  and  simulation  tasks, 
involving  initiation  from  steady-state  operating  conditions  on  the  alternator 
and  loading  circuit  of  various  field,  load,  and  fault  conditions.  For  example, 
partial  comparison  of  percent  deviation  of  analytical  with  respect  to 
experimental  results  for  stabilized  alternator  phase  voltage  and  current 
quantities  after  application  of  transient  load  cases  is  as  follows: 


Load 

Phase  Voltage, 

Phase  Current, 

Case 

percent 

percent 

(a) 

-2.8 

-3.5 

(h) 

-12 

-1.25 

(c) 

-7.2 

1 

(d) 

7 

-12.5 

(e) 

0 

4 

The  above  comparison  data  are  taken  from  Figures  29,  30,  31,  32,  and  3A 
herein  and  load  cases  referred  to  are; 

(a)  Step  increase  of  field  voltage  from  30  percent  to  ratea  with  full 
load  on. 

(b)  Step  decrease  of  field  voltage  from  rated  to  30  percent  with  full 
load  on. 

(c)  Initiate  no  load  to  full  load  at  rated  field  voltage. 

(d)  Single-phase,  1 ine-to-neutra  i sudden  short  circuit  from  rated 
load  at  rated  field  voltage. 

(e)  3"Phase  to  ground,  sudden  short  circuit  from  rated  load  at 
rated  field  voltage. 

Based  upon  the  results  obtained  in  this  work,  it  is  concluded  that  the 
AiResearch  alternator  design  program  ROUND  does  develop  machine  saturation 
data,  resistances,  and  constants  that  may  be  utilized  confidently  to  predict 
alternator  transient  behavior.  Furthermore,  it  is  concluded  that  the  alternator 
dynamic  modeling  and  simulation  mode,  using  either  the  SCEPTRE  or  the 
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AiResearch  computer  programs,  does  give  valid  and  resonably  accurate  results 
compared  with  actual  test. 
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SECTION  II 

DESCRIPTION  OF  REQUIRED  ALTERNATOR  MODIFICATIONS 

The  alternator  selected  for  use  in  this  work  was  a 6-kVA  aircraft  type, 
originally  manufactured  by  Leland  Corporation  and  designated  AGH  329-3.  This 
unit  was  selected  primarily  because  of  its  small  kVA  rating  (which  minimized 
overall  size  and  rating  requirements  for  the  test  circuit)  and  because  it  was 
readily  available.  Modifications  to  the  machine  were  needed  to  make  its 
configuration  similar  to  the  full-scale  machine  designed  under  Contract 
F296OI-75-C-OO55 : however,  a condition  for  using  the  machine  was  that  it  would 
be  restored  to  original  condition  at  the  conclusion  of  the  work. 

The  basic  machine  modification  was  the  use  of  a solid  rotor  oT  nonsalient 
pole  (distributed  winding)  configuration  to  replace  the  original  laminated 
salient  pole  (nond  istr  ibuted  coll)  rotor  structure.  The  shaft  bearing  journals 
of  the  rotor  were  made  wide  enough  to  allow  interference  fit  rings  to  be  pressed 
on  the  bearing  journals  to  permit  loading  the  bearing  races.  The  original  rotor 
had  a center  bore  with  a drive  shaft  extending  the  length  of  the  rotor  shaft  so 
that  the  thrust  bearing  inner  races  could  be  preloaded  by  torque  applied  to 
a bolt  located  at  the  threaded  end  of  the  shaft,  but  with  the  inclusion  of  slip 
rings  on  the  replacement  nonsalient  pole  rotor,  a center-bore  design  could  not 
be  used. 

Slip  rings  and  a brush  setup  were  added  to  allow  the  field  voltage  and 
current  to  be  directly  monitored  under  the  transient  conditions  of  test.  Since 
the  alternator  is  a brushless  design  that  has  a built-in  exciter  machine,  no 
space  was  available  to  incorporate  the  brush-holder  ^issembly  inside  the  stator 
housing  without  major  modification  to  the  stator  assembly.  The  slip  ring/brush 
assembly  was  therefore  located  opposite  the  driving  end  of  the  alternator  and 
was  sized  so  that  the  ring  would  pass  through  the  inner  race  of  the  bearings. 

The  ring  assembly  extended  externally  past  the  end  of  the  stator  housing.  The 
brush-holder  assembly  was  bolted  onto  the  end  of  the  existing  stator  housing. 
Leads  connecting  the  slip  ring  to  the  field  circuit  passed  along  a slot  cut  in 
the  shaft  under  the  bearing  inner  race.  These  modifications  allowed  the 
existing  stator  to  be  used  without  major  structural  or  electrical  changes. 
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The  final  modification  was  the  addition  of  an  air-gap  voltage  search  coil 
to  the  armature.  The  coil  was  wound  in  the  first  slot  of  the  phase  belt  for  one 
phase.  Two  coils  were  wound  in  this  manner  for  two  different  phases.  The 
second  coil  was  redundant  and  was  to  be  used  only  if  the  first  coil  malfunctioned. 
The  coil  is  shown  installed  in  Figure  1. 
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Figure  I.  Air  Gup  Voltage  Search  Coil  Installed  in  Stator 
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SECTION  III 
ALTERNATOR  DESIGN 

Use  of  an  existing  stator  without  winding  or  other  modifications  piaced 
inherent  restrictions  on  the  attainable  reactances  of  the  machine.  (The  only 
modification  made  in  the  stator  element  was  the  addition  of  an  air-gap  voltage 
search  coil.)  Machine  reactances  are  used  to  determine  the  inductances  for  the 
analytical  model  and  many  of  these  are  funct  ions  of  the  stator  geometry  and 
design  (e.g.,  tooth  thickness,  slot  dimensions,  end  turn  length  and  extensions, 
winding  turns,  etc.)  ; however,  basic  changes  do  result  in  machine  subtransient 
and  transient  reactance  quantities  due  to  the  rotor  replacement  made.  With 
careful  consideration  in  design  of  the  field  circuit  and  slot  winding,  the 
solid-core  rotor  without  squirrel  cage  damper  in  the  rotor  slots  does  to  some 
degree  increase  these  reactances. 

The  selection  of  rotor  material  required  a steel  possessing  good  magnetic 
characteristics,  that  is,  capable  of  70  to  90  kilolines  average  rotor  tooth 
flux  densities  without  severe  saturation  and  with  reasonable  structural  proper- 
ties- The  steei  selected  was  4130,  which  exhibits  good  saturation  character- 
istics with  good  mechanical  capability,  is  readily  available,  and  requires  no 
special  material  handling,  tools  or  fixtures,  or  special  heat  treacnent. 

The  original  machine  was  examined  in  detail  to  determine  dimensions, 
stator  winding  parameters,  etc.  AiResearch  computer  program  ROUND  was  then 
utilized  to  develop  the  machine  design  including  optimization  of  the  new 
solid-rotor  conf ' gurat ion . Basic  alternator  performance  data  also  were 
predicted  by  ROUND,  including  the  no-load,  full-load,  and  short-c i rcu I t 
saturation  data  shown  In  Figure  2.  The  major  parameters  for  the  resulting 
modified  design  are  summarized  in  Table  1. 
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Figure  2.  Computer  Result  for  No-Load  and  Full-Loaa  (Constant 
Rated  Current)  Saturation  and  Short-Circuit 


TABLE  I 

MAJOR  PARAMETER  SUMMARY  OF  MODEL  ALTERNATOR 


Stator  Input  Data 


Symbo 1 


Def  i n 1 1 ion 


Q,uant  i 


VBASE 

ABASE 

PHASE 

POLES 

Al 

D 

DSTAK. 

G 

Y 


Base  phase  'citage 
Base  phase  amperes 
Phases 
Poles 

Stator  stack  length,  in. 

Stator  Sore.  in. 

Stator  00,  in.  (either  SETAC  or  DSTAK  must  be 
A i r gap  , In. 

Coil  span,  teeth  or  per  unit 
Total  number  stator  slots 
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TABLE  1 (Continued) 

4 

Symbol 


CHORI 


CVERT 


Stator  Input  Data 


Definition  (Quantity 


Slot  stacking  tolerance,  height,  in.  0 

Slot  stacking  tolerance,  width,  in.  0 

Number  of  parallel  circuits  i .0 

Round  or  rectangular  strand  height  over  insulation,  in,  0.04 

Rectangular  strand  width,  in,  (set  0.0  for  round  wire)  0 

Number  of  horizontal  strands  per  slot  (set  0.0  for  round  0 
w i re) 

Number  of  vertical  strands  per  slot  (set  0.0  for  round  0 

wi  re) 

Number  of  conductors  per  slot  12.0 

Number  of  strands  per  conductor  2.0 

Net  sq  in. /strand  0.00102 

Slot  opening,  in.  0.08 

Tooth  width  (BT  = 0 except  for  straight  tooth)  0.12 

Slot  width,  in.  (Bl  = 0 for  straight  tooth  or  rectangular  0 

wire) 

Height  of  BO , in.  0.02 

Height  of  taper,  in.  0,02 

Height  of  wedge,  in.  0.03 

Thickness  of  slot  liner,  in,  0.015 

Thickness  of  slot  middle  stick,  in.  0.015 

Q.uotient  of  Q./poles /phase  3-0 

S lot  fi  1 1 factor  0.24 

Stacking  factor  0.54 

Current  distortion  factor  (=1.0  for  sine  current)  1.0 

Slot  s kew , I n . 0.31? 

Desired  stator  core  flux  density,  kilolines/sq  in.  0 

Stator  temperature,  °F  350.00 

Alternator  (0.),  Motor  (1.)  0 

Total  length  of  stack  vents,  in.  0 

Densitv,  stator  and  rotor  Iron.  Ib/cu  in.  I 0.276 
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TABLE  I (Continued) 


Stator  Input  Data 

j S '/tnbo  1 

Definition  j 

Q,uant  1 ty 

||BH 

Flux  density,  kl/in^ 

77.0 

Frequency,  Hz 

400.0 

BBli 

Iron  Loss,  Watts/ lb 

10.0 

Rotor  Input  Data 

RSI 

Rotor  slot  grip  height,  in. 

0.06 

RS2 

Rotor  slot  wedge  height,  i i. 

0.02 

RS3 

Rotor  slot  damper  height,  in. 

0. 

RS4 

Rotor  slot  liner  thickness,  in. 

0.02 

RS5 

Rotor  slot  cushion  height,  in. 

0. 

RS6 

Rotor  conductor  insulation  thickness,  in. 

0.001 

RS7 

Rotor  conductor  height,  in. 

0.05 

RS8 

Rotor  bottom  stick  height,  in. 

0 

RS9 

Rotor  conductor  width,  in. 

0.1 

RSW 

Rotor  width , in. 

0.14 

SFCU 

i 

Slope  and  tolerance  for  slot  height,  in. 

0. 

j Q.W 

Number  of  wound  slots  per  pole 

6. 

QS 

Number  of  vent  slots  per  pole 

0. 

Q.TM? 

Number  of  rotor  slots  if  all  slots  were  punched 

8. 

SFPOL 

Rotor  LAM  factor 

1 . 

CONDF 

Number  of  field  conductors  per  slot 

12.00 

CCNCU 

Net  sq  in. /field  conductor 

0.004 

DRI 

Rotor  ID,  in. 

0.0001 

ROTRL 

Rotor  axial  magnetic  length,  in. 

1.85 

EXTRV 

Field  extension  for  venting,  in. 

0.0 

BENDR 

Field  conductor  bend  radius,  in. 

0.25 

TEMPF 

Temperature  of  field  copper 

350.00 

PFWl 

1 

1 

Friction  and  windage  at  FIRON.  watts  (set  = -1 . 
to  calculate  power  loss  due  to  friction) 

-1  .00 

1 FSTRA 

Percent  stray  loss  at  ABASE,  percent 

1 .00 

1 PPOLI 

Pole  head  loss  at  FIRON 

21.00 

- — — ‘ 
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TABLE  1 (Continued) 


Rotoi'  Input  Data  | 

S ymbo  1 

Definition 

Q,uan  City 

Derived 

Constants 

! 

GE 

Equivalent  gap,  in. 

0.01 77 

DR 

Rotor  diameter,  in. 

3-6  i 

81 

P lot  width , in . 

0,2036  1 

HI 

Slot  height  above  wedge,  in. 

0.3776 

HS 

Total  slot  height,  in. 

0.4474  1 

KR 

Carter's  coefficient  for  rotor  slots 

1 .34  ! 

ACOND 

Copper  area  for  a stator  conductor,  sq  in. 

0.0020  ! 

EXT 

Stator  coil  extension,  in. 

1 .64 

ELE 

Length  of  conductor  length 

4.67 

ELC 

Stator  conductor  length 

6.296 

TAUP 

Pole  pitch  at  D,  in. 

2.85  1 

TAUS 

Stator  slot  pitch  at  D,  in. 

0.317 

! TAURO 

1 

Rotor  slot  pitch  at  DR,  in. 

0.353  1 

RSD 

Rotor  slot  depth,  In. 

0.751  ! 

AVCON 

Average  length  of  a field  conductor 

3.64  : 

CM 

Demagnetizing  factor 

1.00  1 

CP 

Field  form  ratio,  average/maximum 

0.528  1 

CF 

Field  se 1 f- 1 i nkage  factor 

C.^35  ! 

Cl 

Fundamental  flux  factor,  peak  fundamenta 1 /peak  total 

0. 866 

SCU 

Stator  slot  fill  factor 

0.24 

TALPF 

Temperature  factor,  field 

1,61  1 

TALPS 

Temperature  factor,  stator 

1 .61 

RA 

Resi stance  per  phase  at  TEMPS,  ohms  * 

0.4883  i 

RF 

Field  resistance,  all  coils  and  poles  in  series 
at  TEMPF,  ohms 

0.2881  1 

i 

1 

TPOO 

Lp/Rp.  = open  circuit  time  constant,  sec 

O.O6S7  ’ 

’ WKP 

Pitch  factor,  stator  winding 

0.985 

I WKD 

1 

Distribution  factor,  stator  winding 

0.959  1 

i WKS 

Skew  factor,  stator  slot 

0.995 

1 WKPOS 

Product  of  WKP,  WKD,  and  WKS 

0.9404  : 

I -/LAME 

Specific  permeance,  stator  end  coi’s 

10.24 
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TABLE  I (Continued) 


Rotor  Input  Data 

Symbol 

Oef  i n i t ion 

'luan  t i cy 

ELF 

Field  inductance,  henries 

0.01  98 

YLAMA 

Specific  permeance,  gap 

327. 

YLAMI 

Specific  permeance,  stator  slot 

3.34 

SLAMR 

Specific  permeance,  rotor  slot 

4.2 

ELAMR 

Specific  permeance,  rotor  end  faces 

1 .1 

RLAFE 

Specific  permeance,  field  end  coils 

5.9 

DLAMD 

Specific  permeance,  damper  slot 

3.8 

RLAME 

Specific  permeance,  damper  end  coils  (double  layer, 
full  cage) 

5-9 

FBS 

Barne's  factor,  stator 

2.4 

FBF 

Barne's  factor,  field 

1.9 

FBD 

Barne's  factor,  damper 

1.54 

FKS 

Same's  factor,  stator 

1 .0 

FKF 

Barne's  factor,  field 

0.7 

FKD 

Barne's  factor,  damper 

0.89 

ABI 

Stator  core  area,  net  sq  i n. 

0.55 

HC 

Stater  core,  height,  In. 

0.5674 

DSTAK 

Stack,  00.  in. 

5.66 

DEN  EM 

WTOT  X 0.785  X DSTAK“  x A! 

0.3 

WCORE 

Stator  core,  weignt,  lb 

3.83 

WFLD 

Field  copper  weight,  lb 

1 .4  ' 

WRO 1 R 

Rotor  i ron  wei ght , 1 b 

4.00 

• 

.vscu 

Stator  copper  weight,  lb 

2.44  ! 

WTETH 

Stator  teeth  weight,  lb 

0.9 

WTOT 

Tota  I electromagnetic  weight,  Ib 

12. 

Load  Condition 

VPU 

AC  volts,  PU 

1 . 

APU 

AC  current,  PU 

1. 

PF 

Power  factor 

-1 . 

RPM 

Revolutions  per  minute 

12000. 

EFF 

Efficiency,  PU 

0.803  i 

Volts 

Phase  voltage 

115.  ; 
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TABLE  1 (Continued) 
Load  Cond i t ion 


I 


S ymbo 1 


Def in i t ion 


Q,uant  i tv 


AMPS 

Phase  current 

17.^ 

KVA 

KVA  rati ng 

6. 

FREQ 

Frequency,  Hz 

400. 

KW 

Load , kW 



Variable  Derived  Parameters 

AMPSF 

Fi eld  current,  A dc 

29.2 

CHI 

Reactance  factor 

0.0137 

DENFD 

Current  density,  field 

7301  . 

DEN  ST 

Current  densi  ty,  stator 

8529. 

ELOAD 

Electrical  loading,  A-conductor/i n. 

619. 

HEATF 

ELOAD  X DEN  ST 

5.3  X 10^ 

PHI  LK 

Field  leakage  flux  per  pole,  kilolines 

8.7 

PHI  PL 

Air  gap  flux  per  pole,  kilolines 

TOO. 

aUALF 

kVA/lb/ rpm 

4.5  X 10“^ 

RAPU 

RA/Z  base,  per  uni t ohms  (hot) 

0.074 

< 

1- 

Armature  time  constant  for  the  dc  component, 
X^/Z  FRAPU,  sec 

0.001 2 

THETA 

Power  factor  angle,  radian 

0. 

TPO 

Transient  time  constant,  sec 

U^ 

o 

O 

o 

Z3ASE 

Base  resi stance,  ohms 

6.6 

BROOT 

Flux  density  at  root  of  rotor  tooth  next  to  pole 
center,  kilolines/sc  in. 

127. 

3ETAC 

Flux  density,  stator  core,  kilolines/sq  in. 

43. 

BETAG 

Flux  density,  air  gap,  kilolines/sq  in. 

40.7 

. 

BPCNT 

Flux  density,  pole  center,  kilolines  sq  in. 

46.7 

BETAT 

Flux  density,  stator  teeth,  kilolines/sq  in. 

114.00 

1 

BETCO 

Flux  density,  rotor  core,  kilolines/sq  in. 

28. 

FCORE 

Magnetomotive  force,  stator  core,  A-turns 

3.4 

J 

FGAP 

Magnetomotive  force,  air  gap,  A-turns 

225. 

FPCNT 

Magnetomot i ve  force,  pole  center,  ^-tur-'S 

29. 

FTETH 

Magnetomotive  force,  stator  teeth,  A-turns 

142, 

1 

FRO  CO 

MacnetonxDt ! ve  force,  ro^or  core.  A-turns 

■ 

1 A. 

1 

15 


AFWL-TR- 75-66,  Add.  2 


r 


TABLE  1 (Continued) 


Variable  Derived  Parameters 


S ymbo 1 

Definition 

Quan  c i t V 

FIELD 

Magnetomotive  force,  field  per  pole,  A-turns 

1051.  j 

FFL 

Magnetomotive  force,  behind  XP , A-turns 

862. 

HRC 

Height,  rotor  core,  in. 

1 .04 

3RC 

2 

Flux  density,  maximum  in  rotor  core,  kilolines  in. 

31  . 

PCORE 

Stator  core  loss,  watts 

31. 

PFLDT 

Field  copper  loss,  watts 

245. 

PFV 

Friction  and  windage  loss,  watts 

39. 

PPOLH 

Pole  head  loss,  watts 

665. 

PSCU 

Stator  copper  loss  , watts 

443. 

PSTRA 

Stray  load  loss,  watts 

3.5 

PTtTH 

Stator  teeth  loss,  watts 

47. 

PTOT 

Total  loss,  watts 

1476. 

XAO 

Direct  armature  reaction,  PU 

3.84 

X2 

Negative  sequence  reactance,  PU 

0.219 

XD 

Direct  unsaturated  synchronous  reactance,  PU 

4.028 

XDDL 

Direct  damper  leakage  reactance,  PU 

0.0266 

XO 

Zero  sequence  reactance,  PU 

0.085 

XF 

Field  leakage  reactance,  PU 

0.195 

XL 

Stator  leakage  reactance,  PU 

0.184 

XPO 

Direct  transient  reactance,  saturated,  PU 

0.323 

XPDU 

Direct  transient  reactance,  unsaturated,  ?U 

0.367 

XPPD 

Direct  subLransient  reactance,  PU 

0.21  0 

XDQ 

Quadrature  damper  reactance,  PU 

3.3S 

XDQL 

Damper  quadrature  leakagt^,  reactance,  PU 

0.0449 

XPPQ. 

Quadrature  subtransient  reactance,  PU 

0.228 

TPPD 

Direct  subtransient  time  constant,  sec 

0.00525 

XAQ. 

Quadrature  armature  reaction,  PU 

3.31 

Xd 

Quadrature  synchronous  reactance,  PU 

3.49 

XDD 

Direct  damoer  reactance,  PU 

3.87 

XPO. 

Quadrature  transient  reactance,  PU 

0.2286 

1 


A 
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SECTION  IV 

ALTERNATOR  FABRICATION,  ASSEMBLY,  AND  TEST  SETUP 

Figure  3 shows  the  replacement  rotor  after  milling  slots.  The  smallest 
diameter  shaft  extention  was  later  knurled  to  retain  the  slip  ring  assembly. 

The  slotted  diameter  is  the  bearing  journal  through  which  the  lead  connecting 
the  slip  rings  to  the  field  circuit  is  passed.  On  the  opposite  end  of  the 
core  the  long  shaft  extension  carries  the  drive  end  bearing  and  drive  coupling 
provisions.  This  area  of  the  shaft  in  the  original  brushless  (salient  pole) 
rotor  held  the  rotating  rectifier  assembly  armature  for  the  main  exciter. 

The  final  assembled  rotor  is  shown  in  Figure  4.  The  slip  ring  assembly 
is  shown  fitted  to  the  rotor.  The  slip  rings  were  separated  and  insulated 
from  the  rotor  shaft  by  Astrel  360  (pol yaryl sul fone)  plastic,  which  is  a 
good  dielectric  that  provides  mechanical  stability  and  strength  and  high- 
temperature  capability. 

A 0. 01 0- in. -thick  slot  liner  of  Nomex-Kapton-Nomex  (NKN)  insulator  was  | 

used  to  insulate  the  pole  turns  from  the  rotor  core.  Each  turn  consisted  of  | 

two  parallel  No.  16  AWG  round  magnet  wires  having  polyamide-imide  insulation  | 

meeting  Class  220  (220^C  continuous)  thermal  rating.  To  achieve  the  maximum 
slot  fill  factor,  the  turns  of  two  No.  16  wires  were  laid  side  by  side  in  the  J 

slot  without  crossovers.  As  each  coil  was  completed  the  slot  liner  was  folded 
over  the  turns  and  310  nonmagnetic  stainless  steel  top  wedges  were  inserted  to 
retain  the  coils  in  the  slot. 

After  all  poles  had  been  wound,  the  coils  on  each  pole  were  series-  ; 

connected.  All  pole  groups  were  then  series-connected  and  solder-terminated  to 
the  field  leads.  The  end  turns  on  each  pole  were  banded  together  using  a 
prestressed  fiberglass  tape,  and  a prestressed  band  of  fiberglass  tape  was 
then  wound  around  the  OD  of  all  pole  end  turns  to  restrain  them  against 
centrifugal  loading  forces.  The  fiberglass  tape  was  cured  for  4 hr  at  325°F. 

The  slip  ring  assembly  was  then  pressed  onto  the  rotor  shaft  and  the  adapter 
for  the  drive  coupling  was  copper  flash- coated  and  press-fitted  into  the  drive 
end  of  the  rotor.  All  surfaces  were  then  masked,  and  the  rotor  was  vacuum 
impregnated  in  silicon  varnish.  The  cure  cycle  for  the  varnish  was  5 hr  at 
325°F. 

J 
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Insulation  resistance  tests  were  conducted  at  500  vdc  during  the  winding 
procedure  as  each  coil  group  was  completed  and  again  at  end  of  each  cure  cycle. 
No  breakdown  in  rotor  conductor  insulators  occurred.  The  dc  resistance  of  the 
rotor  winding  was  measured  at  completion  and  was  0.195  ohm  at  70°F. 

The  completed  rotor  assembly  was  balanced  including  bearings.  Balance 
was  accomplished  by  drilling  and  grinding  away  rotor  core  and  shaft  material. 

A final  unbalance  tolerance  of  0.003  in.-oz  was  achieved. 

The  rotor  was  placed  in  the  stator  after  installing  the  stator  search 
coils.  The  two  end  bearing  race  loading  rings  were  interference-fitted  on  the 
shaft,  and  bearing  inner  races  were  loaded  by  pressing  from  both  ends  of  the 
alternator.  The  machine  was  then  reassembled  and  the  brush-holder  assembly 
containing  four  electrographi tic  brushes  of  0.5  in.  by  0.25  in.  surface  area 
was  bolted  to  the  stator  housing.  Brushes  were  contoured  to  the  slip  ring 
surface  by  means  of  fine  sand  paper  placed  on  the  slip  ring  while  slowly  rota- 
ting the  shaft.  Care  was  taken  to  remove  any  remaining  abrasive  material  from 
the  slip  ring  and  brush  holder.  The  assembled  structure  was  then  coupled  to 
the  electrical  dynamometer  as  shown  in  Figure  5. 

The  alternator  was  driven  by  the  electrical  dynamometer,  which  was  struc- 
turally connected  to  i t by  means  of  the  alternator  end  bell  flange  bolts.  The 
power  output  of  the  dynamometer  is  9 hp  at  12,000  rpm.  The  alternator  shaft 
was  connected  to  the  dynamometer  drive  shaft  through  a specially  constructed 
flexible  shaft  coupling.  The  alternator  is  a forced-air-cooled  design,  and 
cooling  air  was  provided  by  a test  facility  air  system.  The  cooling  air  was 
introduced  through  the  existing  ducts.  Stator  output  power  and  field  leads 
(from  the  brush  assembly)  were  connected  to  Jones  strip  connectors  located  on 
a metal  plate  bolted  to  the  alternator  mounting  base.  This  allowed  easy 
removal  and  change  of  test  leads  during  the  test.  Temperature  sensors  and  the 
voltage  search  coils  incorporated  into  the  stator  housing  also  were  terminated 
on  the  Jones  strip  connector. 
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SECTION  V 

TEST  CIRCUIT  DESCRIPTION 

Figure  6 is  a schematic  of  the  alternator  loading  circuit  constructed. 

Each  of  the  three  phases  of  the  alternator  was  connected  to  four  capacitors 
(12  total)  arranged  so  that  the  two  capacitors  on  the  positive  half  cycle  and 
two  capacitors  on  the  negative  half-cycle  of  each  phase  were  charged.  The 
charging  of  the  capacitors  was  sequential  and  was  controlled  by  turning  the 
SCR's  on  or  off  at  proper  intervals.  The  gate  firing  logic  circuit  for  each 
SCR  determined  its  on  or  off  state.  There  were  12  charge  SCR's  ((11C-Q12C)  and 
12  discharge  SCR's  (Q10-Q12D). 

The  gate  firing  logic  was  designed  so  that  each  charging  SCR  was  turned 
on  once  in  two  full  cycles,  and  the  SCR  stayed  on  for  one-half  cycle.  In  a 
similar  manner,  each  discharge  SCR  was  turned  on  once  in  two  cycles,  and  stayed 
on  for  approximately  one- th i rd  cycle. 

Figure  7 shows  the  capacitor  charge  and  discharge  currents  at  ideal 
(resonance)  conditions.  The  resonant  charging  frequency  is  determined  by  the 
capacitance  and  effective  phase  inductance  value  of  the  alternator.  The  dis- 
charge period  of  the  capacitor  was  determined  by  the  RC  time  constant  of  the 

capacitor  and  load  resistor  adjusted  to  allow  five  RC  time  constants  for  the 

total  discharge  of  the  capacitor  in  the  one-third  cycle  time  period. 

The  actual  test  parameter  component  values  for  the  resonant  load  capaci- 
tors and  the  load  resistors  are  listed  in  Table  2. 

Figure  8 is  a block  diagram  for  the  load  firing  control  circuit  of  Figure  5 

The  following  functional  description  explains  the  control  from  tiie  input  of  the 

air  gap  search  coil  voltage  to  the  outputted  sequential  firing  of  the  charge 
and  discharge  SCR's. 

The  Frequency  Reference  Generator — The  alternator  search  coil  provides 
a sync  signal,  that  is,  the  air  gap  voltage  to  the  frequency  refer- 
ence generator.  The  sync  signal  is  then  phase- d i sp 1 aced  130  deg  to 
eliminate  noise.  The  phase-displaced  signal  generates  a square  wave 
signal  whose  frequency  is  the  same  as  the  frequency  of  the  search  coil 

Phase  Lock  and  Clock  Generator — The  scuare  wave  output  of  the 
frequency  reference  generator  is  compared  with  a frequency 
feedback  signal  from  the  phase  generator.  This  provides  phase 
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locking  between  the  two  signals  and  results  in  a clock  pulse  to  be 
generated  whose  frequency  is  six  times  the  square  wave's  frequency 
from  the  frequency  reference  generator. 

The  Phase  Generators— The  clock  pulse  is  inputted  to  three  flip-flops 
connected  to  provide  three  phase  logic  signals  called  A,  8,  and  C, 
each  120  deg  apart  and  their  inverses  called  A,  B,  and  C.  The  A,  3, 
and  C logic  signals  are  fed  back  to  complete  the  phase  lock  control 
1 oop. 

Gate  Driver  Logic— The  three  phase  logic  signals  A,  B,  and  C and  the 
three  inverted  logic  signals  A,  B and  C from  the  two  phase  genera- 
tors are  terminated  into  the  gate  driver  logic  to  generate  firing 
command  signals  to  the  SCR,  Each  phase,  that  is  A and  A,  provides 
the  driver  logic  signal  command  to  charge  and  discharge  sequentially 
four  capacitors  according  to  the  order  shown  in  Figure  7. 

Field  Gate  Drivers--The  driver  logic  signals  from  the  gate  driver 
logic  are  terminated  by  eight  field  gate  driver  circuit  cards. 

The  drivers  provide  the  necessary  gate  driver  power  signal  to 
fire  the  SCR's.  Each  circuit  card  contains  drivers  for  three  SCR 
gates. 

Figure  9 shows  the  actual  test  circuit  layout.  The  SCR  heat  sink  assembly 
and  the  logic  card  rack  rest  upon  the  top  of  the  table  and  the  capacitor  bank 
underneath. 
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SECTION  VI 
TEST  DESCRIPTION 


The  test  alternator  program  encompassed  four  successive  tasks; 

• Basic  machine  magnetic  tests 

• Determination  of  alternator  resistances  and  reactances 

• Tests  with  the  load  circuit  under  steady-state  conditions 

• Tests  with  the  load  circuit  under  transient  conditions 

Tests  conducted  to  determine  alternator  reactances  and  saturation  data  were 
performed  in  accordance  with  the  IEEE  No.  115  Test  Procedures  for  Synchronous 
Machines.  Specifically  these  were: 


Test 

Open-circuit  saturation 
Short-circuit  curve 

Di rsct-axi s transient  and  subtransient 
Reactance  test  from  sudden  short-circuit 
Direct  synchronous  reactance 
Direct  and  quadrature  subtransient  reactance 
Zero-sequence  resistance  and  reactance 


Test  Section 

3.05.35 

3.05.50 

7.25.15  Method  1 
7.30,10  Method  1 

7.15 

7. 30. 25  Method  3 
7.  ^5.  1 0 Method  1 


TEST  EQUIPMENT  AND  MONITORED  PARAMETERS 

The  test  equipment  and  test  parameters  are  summarized  in  Table  3.  Figures 
10  and  11  show  the  test  equipment  physical  layout  including  power  supolies, 
transient  recorder,  and  metering  equipment  used  in  the  testing  with  the  loading 
circuit. 

BASIC  MAGNETIC  CHARACTERISTICS 

Basic  machine  characteristics  were  obtained  at  (1)  open-circuit,  (2)  short- 
circuit,  and  (3)  full-load  conditions.  The  full-load  tests  were  performed  using 
a constant  value  3-phase  resistive  load.  The  full-load  saturation  curve  was 
determined  in  a manner  similar  to  that  for  the  circuit  saturation  by  terminating 
the  alternator  with  a 3-phase  resistive  load,  which  at  rated  full  line  to  neutral 
voltage,  develops  full-load  phase  current  and  power  output.  Test  results  are 
shov/n  in  Figure  12. 
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TEST  OF  ALTERNATOR  RESISTANCES  AND  REACTANCES 

Test  of  alternator  parameters  included  the  following  tests: 

(a)  Direct-axis  transient  and  subtransient  reactance  from  sudden  short- 
c i rcu  i t. 

(b)  Direct  and  quadrature  subtransient  reactance 

(c)  Zero  sequence  reactance  and  resistance 

(d)  Direct  synchronous  reactance 

The  following  summarizes  the  per  unit  (PU)  reactance  values  determined 
experimentally  in  the  test: 


Test 

Xd' 

Xd" 

Xq" 

Ro 

Xo 

Xd 

1 

0.449 

0.  201 

— 

— 

— 

— 

2 

— 

0.  21  75 

0.2451 

— 

— 

— 

3 

— 

— 

— 

0.0451 

0.0958 

— 

k 

— 

- — 

— 

— 

— 

2.  72 

^base  = ohms 


TESTS  WITH  THE  LOAD  CIRCUIT  UNDER  STEADY-STATE  CONDITIONS 

The  steady-state  operation  load  circuit  tests  were  conducted  prior  to  the 
initiation  of  the  transient  load  case  tests  in  order  to  determine  if  the  search 
coil  air  gap  voltage  was  sufficient  to  properly  operate  the  load  firing  control 
circuit.  It  was  found  that  the  air  gap  voltage  used  for  triggering  at  rated 
field  voltage  and  speed  achieved  proper  initiation  of  the  charge  and  discharge 
of  capacitors.  The  results  of  the  steady-state  test  are  shown  in  Figures  13 
through  16. 

TESTS  WITH  THE  LOAD  CIRCUIT  UNDER  TRANS  I ENT-COND 1 T 1 ONS 

Table  k summarizes  the  transient  load  case  performance  tests.  In  load 
cases  1 and  2,  field  voltage  could  not  be  stepped  from  zero  to  rated  because 
sufficient  air  gap  voltage  could  not  be  generated  by  the  excited  field  to 
develop  the  sync  lock  feedback  previously  described  in  the  test  circuit  descrip- 
tions. It  was  found  that  a minimum  of  approximately  30  percent  of  full  ratea 
field  voltage  was  necessary  to  maintain  the  sync  lock  under  a step  change  in 
field  vol tage. 
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CALIBRATION:  PHASE  CURRENT : 20  A/CM 

PHASE  VOLTAGE:  50  V/CM 

TIME : 0 5 MS/CM 


Figure  13.  Alternator  Phase-Voltage  and  Current 


CALIBRATION:  PHASE  CURRENT:  20  A/CM 

SEARCH  COIL  VOLTAGE:  5 V/CM 

TIME  : 1 .0  MS/CM 


Figure  14,  Alternator  Phase  Current  and  Search  Coil 


r-2371^ 


Vo  1 tage 
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CALIBRATION:  FIELD  CURRENT:  6 A/CM 

FIELD  VOLTAGE:  2 V/CM 

TIME:  1.0  MS/CM 


ZERO 

LEVEL 


Figure  ^S•  Field  Voltage  and  Current 


CALIBRATION:  LOAD  CURRENT:  20  A/CM 

LOAD  VOLTAGE:  20  V/CM 

TIME:  0.5  MS/CM 


PEAK 

CURRENT 


PEAK 
vni  TAGE 

r-?3/17 


Figure  16.  Load  Voltage  and  Current 
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TABLE  k 

LOAD  CASES  FOR  ALTERNATOR 


Case 

No, 

i 

Torque 

Speed 

Field  Voltage 

Load 

1 

Rated 

Rated 

Step  increase  from 
30  percent  to  rated 

1 

Rated  1 

2 

Rated 

Rated 

Step  decrease  from 
rated  to  30  percent 

Rated 

1 

i 

3 

Rated 

Rated 

Rated 

1 

Step  increase  from  zero  to  ; 

rated  i 

A 

Rated 

Rated 

Rated 

1 -phase,  1 i ne-to-neut ra 1 short! 
circuit  from  rated  load  j 

5 

Rated 

Rated 

Rated 

1-phase,  1 i ne- to- 1 i ne,  short  j 
circuit  from  rated  load  i 

6 

Rated 

Rated 

Rated 

3-phase  to  ground,  short 
circuit  from  rated  load  ^ 

7 

Rated 

Rated 

Rated 

Open-circuit  from  rated  load,  1 
1 -phase  1 i ne-to-neutra 1 | 

The  test  results  are  presented  in  Figures  17  through  23.  The  calibration 
of  the  recordings  is  as  follows: 


Gal vanometer 

Number Parameter  Ca1  i brati  on 


1 

A voltage,  1 i ne- to-neutra 1 

260  V/in. 

4 

A phase  current 

30  A/ in. 

11 

B phase  current 

80  A/ in. 

14 

Search  coil  vol tage 

7 V/in, 

15 

Field  current 

10  A/ in. 

18 

Field  voltage 

5 V/in. 

Time:  0.0156  s/in. 
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Lodd  Case  1 Transient  Performance  Test  Results  r-a/u-j 


r ~~ 
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Fiijure  18.  Load  Case  2 Transient  Performar^e  Test  Results 


Load  Case  3 Transient  Performance  Test  Results  f-?37s/ 


ure  20.  Load  Case  4 Transient  Performance  Tes 


Figure  22.  Load  Case  6 Transient  Performance  Test  Results  F-23/iS 
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SECTION  V I I 

ALTERNATOR  MODELING  APPROACH 

The  alternator  model  used  in  this  analysis  is  shown  in  Figure  2A. 


Figure  24.  Alternator  Moae 1 

This  model  consists  of  six  circuits:  three  armature  phase  windings;  one 

field  winding;  one  direct  axis  damper  winding;  and  one  quadrature  axis  camper 
winding.  The  equivalent  damper  windings  usually  have  the  same  number  of  turns 
as  the  armature  wincings.  The  field  winding,  on  the  other  hand,  usually  retains 
its  own  number  of  turns.  There  will  be  mutual  incuctances  between  these  vvinc- 
ings.  The  speed  of  the  machine  is  assumed  constant. 

The  matrix  equation  for  the  six-circuit  model  is  given  by; 


(Ml)  m R|  - V = E (1) 

Since  inductances  are  affected  by  saturation,  Ecua  t i on  (1)  can  be  written  as 


AFWL-TR-75-66,  Add.  2 


r 


Before  discussing  the  solution  of  this  set  of  differential  equations,  it  is 
necessary  to  discuss  the  alternator  constants.  The  AiResearch  alternator 
design  computer  program  ROUND  calculates  the  unsaturated  parameters  (based  on 
an  air  gap  line)  x,.  X^,.  X'^,  X^",  X,^,  X„^,  X„^,,  X„^,.  X^,  X^,  X^, 

and  L^-,  as  well  as  a predicted  no-load  saturation  curve.  From  these  values, 

together  with  rated  frequency  and  the  machine  base  impedance,  the  following 
unsaturated  machine  inductances  can  be  obtained; 

Lj  Direct  axis  inductance,  H 

Lq  ^ Quadrature  axis  inductance,  H 

Armature  leakage  inductance  per  phase,  H 
Direct  axis  transient  inductance , H 
Direct  axis  subtransient  inductance,  H 
Quadrature  axis  subtransient  inductance,  H 
Direct  axis  damper  winding  inductance,  H 
Quadrature  axis  damper  winding  inductance,  K 
Direct  axis  damper  winding  leakage  inductance,  H 
Quadrature  axis  damper  winding  leakage  inductance,  h 
Negative  sequence  inductance,  h 
Zero  sequence  inductance,  H 

Field  winding  leakage  inductance,  in  terms  of 
armature,  H 

Field  winding  se 1 f- i nductance  (in  its  own  term)^  H 


al 

L'd 

L". 


■-Dd 

Sq 

'-Ddl 

‘■Dql 

*■2 

I 

"o 

^'fl 

■-ff 


The  inductances  and  resistance  of  the  alternator  model  are  as  follows. 
(These  inductances  can  be  derived  from  the  direct  and  quadrature  axes  induc- 
tances listed  above). 


'aa 

Self 

i nductance 

of 

a rrr.a  t u re 

w ind  inc 

phase 

A, 

M 

“■bb 

Se  1 f 

i nductance 

of 

armature 

winding 

phase 

3. 

J 

^cc 

Self 

i nductance 

of 

arr.a  ture 

w i nd  i ng 

ohase 

H 
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Mutual  inouctance  between  armature  windings  A and  8,  H 

Mutual  inductance  between  armature  windings  B and  C,  H 

Mutual  inductance  between  armature  windings  C and  A,  H 

Self  inductance  of  field  winding,  in  its  own  term,  H 

Self  inductance  of  direct  damper  winding,  in  armature 
terms,  H 

Self  inductance  of  quadrature  damper  winding,  in 
armature  terms,  H 

Mutual  inductance  between  field  and  armature  winding 
Phase  A , H 

Mutual  inductance  between  field  and  armature  winding 
Phase  B,  H 

Mutual  inductance  between  field  and  armature  winding 
Phase  C,  H 

Mutual  inductance  between  direct  damper  and  armature 
winding  Phase  A , H 

Mutual  inductance  between  direct  damper  and  armature 
winding  Phase  B,  H 

Mutual  inductance  between  direct  damper  and  armature 
winding  Phase  C,  H 

Mutual  inductance  between  quadrature  damper  and 
armature  -winding  Phase  A,  H 

Mutual  inductance  between  quadrature  damper  and 
armature  winding  Phase  8,  H 

Mutual  inductance  between  quadrature  damper  and 
armature  winding  Phase  C,  H 

Mutual  inductance  between  field  and  direct  damper 
w i nd i ngs  , H 

Armature  winding  leakage  inductance,  same  for  all 
three  phases,  H (portion  of  the  armature  inductance 
not  subject  to  saturation) 

Field  winding  leakage  i nductance  ( i n its  own  term),  H 
Direct  damper  winding  leakage  inductance,  H 
Quadrature  damper  winding  leakage  inductance,  H 
Armature  -winding  resistance  per  phase,  ohms 
Field  -winding  resistance,  or, ms 


-sst 


hS 
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Direct  damper  winding  resistance,  ohms 
Quadrature  damper  winding  resistance,  ohms 


The  mutual  inductances  between  the  field  and  the  quadrature  damper  as 
well  as  between  the  direct  and  the  quadrature  damper  windings  are  assumed  to 
be  zero  even  though  the  saturation  effect  would  cause  some  slight  coupling 
between  them. 

The  alternator  under  consideration  has  a round,  solid  rotor.  Uneven 
slotting  gives  a slight  sal iency  to  the  machine.  Damping  effect  comes  from 
eddy  currents  induced  in  the  solid  rotor  structure. 

The  self  and  mutual  inductances  of  the  armature  windings  can  be  approxi- 
mated by  the  following  expressions  (the  sal iency  is  approximated  by  a double- 
frequency sinusoidal  term.) 

L = L + L cos  2 aut  (12) 

aa  s m 


L,,,  = L L cos  (2  -i-t  + 120) 
bb  s m 


L = L + L cos  (2  i;t  - 120) 
cc  s m ' 


M = -M  - L cos  (2  ^t  -t  60) 

ab  ba  s m 


M = -M 
cb  s 


L cos  2 JJt 
m 


M - L cos  (2  JJt  - 60) 
s m 


ca  ac  s m 

where  xt  is  expressed  in  electrical  degrees. 


The  mutual  inductances  between  field  or  damper  windings  with  armature 
windings  vary  with  relative  winding  positions.  This  variation  is  approximated 
by  a sinusoidal  function.  The  positive  direction  of  the  quadrature  axis  is 
taken  as  90°  ahead  of  that  of  the  direct  axis. 


"af  = 

M,  = 
fa 

"af 

cos 

■Ut 

(18) 

^f  = 

rb 

"af 

cos 

(J-t  - 

120) 

■19) 

\f  = 

O 

II 

"af 

cos 

(OJt  -r 

120) 

■20) 

^aOd 

" ^Dda 

^ ^add 

cos  J/t 

^21) 

'"’bOd 

* ^Ddb 

^ ^ADD 

cos  (-tt 

- 120) 

(22) 

47 


AFWL-TR- 75-66,  Add.  2 


^Dd  = '^Ddc  = ^DD  ^ ’20) 


’'aDq  = ’^Dqa  = ’^ADQ  ^ ^O) 

^Dq  = %qb  " "AOd  ‘ 3°) 

^Dq  = ^Dqc  = ^ADd  " 210) 


(23) 

[2k) 

(25) 

(26) 


Other  mutual  inductances  are 


^fOd  = ^Ddf  = ^ <:onstant 


^fDq  " '^Dqf  " ° 


(27) 

^28) 


M = M =0 
DdDq  DqDd 


(29) 


Those  constants  not  listed  in  the  computer  output  are  calculated  as 
fo 1 1 ows ; 


= 1/3  (L,  - - 

^ = 1/2  (L^  - L^) 


’•m  = 2/3  (L,  - - M^) 


AF 


M..  = \/-|  L„  (L. 


■ff  ■‘‘d 


d) 


(L,,  Is  in  field  terms,  not  in  terms  of  the  armature) 

I 1 

Let  k - = effective  turns  ratio,  armiature  (one  phase  only,  not 
d T 

three-phase  equivalent)  to  field. 


30) 

31) 

32) 

33) 


'ADD 


'^af’^AF 


^ADQ  Y 3 ‘■Dq^S"‘‘q^ 

The  value  of  can  be  ootained  by  solving  the  following  quadratic 


(3^) 

(35) 


equat  ion. 


"FDD 


ADD  AF 

tT 

Ld  - Ld 


^DD  ^ 


3 ‘‘Dd'’''AF  ‘■ff''’ftDD  , , 

3 (Tt;^ 


= 0 


(36) 


U8 
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The  smaller  root  is  used,  which  is  the  physical  solution.  This 

2 


quadratic  equation  comes  from  the  relation: 

^ 1 ‘-Pd'^AF  " ^^FDD‘‘^ADD^AF  ‘‘ff^ADD 


L.  - L' 


(37) 


‘-Dd'-ff  " ^FDD 


In  a solid  rotor  alternator  without  pole  face  type  damper  winding,  the 
subtransient  component  decays  very  rapidly.  Using  one  damping  winding  on 
each  direct  and  quadrature  axis,  a good  compromised  open-circuit  subtransient 
time  constant  is  1/3  cycle.  In  a Aoo-Hz  machine. 


and 


T , = T a:  0.00085  sec 

do  qo 


Dd 


■Dd 


do 


Dq 


(38) 

(39) 


qo 

In  this  experimental  alternator  with  round  rotor  construction  it  is  consi- 
dered reasonable  to  use  one  saturation  factor  for  all  inductances.  The 
saturation  factor  is  determined  by  the  resultant  MMF  (field  ampere-turns  plus 
the  armature  reaction). 

The  resultant  MMF  can  be  resolved  into  two  components,  MMF^  along  the  direct 
axis  and  ^MF^  along  the  quadratui'e  axis.  Armature,  field,  and  direct  damper 
currents  contribute  to  MMF^.  Armature  and  quadrature  damper  currents  contribute 
to  MMFq.  The  MMF  producedby  any  armature  phase  current  is  pulsating  in 
magnituae  and  stationary  in  space.  This  MMF  can  be  resolved  into  two  opposite- 
revolving  MMF's,  each  with  one  half  the  magnitude  of  the  pulsating  MMF.  The 
forward  revolving  component,  which  is  stationary  with  respect  to  the  rotating 
field,  is  the  armature  reaction  component.  The  three  backward  revolving 
components  of  the  three  phases  will  cancel  each  other  in  the  case  of  a 
balanced  load.  In  case  of  unbalanced  load,  the  resultant  MMF  will  revolve 
with  double  speed  relative  to  the  field  structure.  The  contribution  of  this 
double-speed  MMF  to  saturat  ion  of  inductances  is  ignored  in  the  analysis. 


i,9 
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If  the  MMF's  are  expressed  in  terms  of  equivalent  field  amperes,  and 
at  time  zero,  the  field  axis  coincides  with  the  axis  of  armature  winding 
phase  A,  the  instantaneous  values  of  MMF^  and  MMF^  can  be  expressed  in  terms 
of  the  instantaneous  currents  as  follows; 

'^''d  = ‘f  ^ '"af’Dq  ^ I ^f  I'a  'b  - 1 20  ) x 

i^  cos  (xt  + 120°)]  (40) 

MMF  = k ,i  „ - ^ k ,fi  sin  xt  + i,  sin  (xt  - 1 20  ) + 

q af  Dq  2 af[  a b 

i ^ sin  (xt  + 1 20°)|  (4l ) 

As  mentioned  previously,  the  positive  direction  of  this  quadrature 
axis  is  assumed  to  be  90  deg  ahead  of  that  of  the  direct  axis.  The  resultant 
MMF  can  then  be  obtained; 


MMF  = /MMF 


MMF 


(42) 


With  unbalanced  load  and  during  transient,  these  MMF's  vary  with 
time.  From  the  no-load  saturation  curve  and  the  air  gap  line,  a saturation 
factor  can  be  defined  from  Figure  25  such  that  at  a certain  resultant  MMF  of 

‘fl. 

F 

sat  AC 

Figure  26  shows  the  result  of  this  calculation. 

Since  in  general,  MMF  is  a function  of  time,  the  saturation  factor  is  also 
a function  of  time. 

The  saturated  inductances  are  obtained  by  multiplying  the  unsaturatea 
i nductance  by  the  saturation  factor.  For  instance,  the  saturated  se 1 f-l nductance 
of  armature  winding  phase  A is; 


aa,sat  al 


*^sat  ^‘■aa  ‘ '"aP 


(44) 


where  is  leakage  inductance  and  is  for  practical  purposes  not  subject 
to  saturation. 

For  the  experimental  alternator,  the  constants  from  the  Ai Research  com- 
puter design  program  ROUND  are  tabulated  in  Tables  5 througn  7- 
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TABLE  5 

COMPUTER  OUTPUT  CONSTANT 


><d 

4.028  PU  1 

X 

3.498  PU  1 

q 1 

^al 

0.1837  PU  : 

^'d 

0.3676  PU  j 

X" 

i ° 

0.2103  PU  1 

X" 

0.2286  PU  i 

q 1 

Sd 

3.S71  PU  1 

Sq 

X.  .. 

3.359  PU 

0.0266  PU 

Ddl  1 

^Dql 

0.0449  PU 

^2 

0.2195  PU 

*0 

0.08535  PU 

^fl 

0.1931  PU 

L^^  = 0.01981 

H 

TABLE  6 


ALTERNATOR  MODEL  CONSTANTS  DERIVED  FROM  COMPUTER  OUTPUT 


L 

0.00667  H 

s 

M 

0.00322  H 

s 

! 

0.000467  H 

i 

O.OI9SI  H 

' L.  . 

0.01018  H 

i Dd 

‘-Dq 

M 

0.00883  H 

0.01128  H 

«F  i 

^ADD 

0.00825  H 

M ... 

0.00712  H 

1 ADli 

M 

0.0142  H ! 

FDD 

1 L , 

C.  000483  H 

ai 

1 '■fl 

0.000948  H 

! Sdi 

0.000070  H 

i Sql 

0.0001 1 8 H 
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TABLE  6 (Conti nued) 


R 

a 

0.300  ohm 

~ " ' 1 

'^f 

0.2665  ohm 

(i nc ludi ng  slip  ri ng 
brush  drop; 

Sd 

12.0  ohm 

’^Dq 

10.0  ohm 

TABLE  7 

SATURATION  FACTOR  DATA  (FOR  COMPUTER  PROGRAM  ROUND  CONSTANTS) 


where 

A = 0.02473334 
B = 7-578884 
C = - 18.77908 


Using  test  results  on  the  experimental  alternator  combined  with  manually 
calculated  constants  using  Kilgore's  design  formulas  (Reference  3),  the 
circuit  constants  for  the  model  are  those  in  Table  8 and  9.  These  constants 
are  used  together  with  a different  set  of  saturation  curves,  from  actual 
test.  As  shown  in  Figure  27,  a fictious  air  gap  line  is  drawn  tangential 
to  the  lower  part  of  the  tested  no*load  saturation  curve.  The  saturation 
factor  at  a MMF  of  1^]  is  again  given  by 


F = ^ 
sat  AC 


Figure  28  shows  the  result  of  Equati on  C45  ) . 


Reference  3-  KILGORE,  L.  A.,  Calculation  of  Synchronous  Macn i ne  Constants-- 

Reactances  and  Time  Constants  Affecting  Transient  Characteristics, 
pp  1201-1214,  AIEE  Transact  ions , December  1931 
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TABLE  8 

ALTERNATOR  MODEL  CONSTANTS  DERIVED  FROM  COMBINATION 


OF  TEST  AND  COMPUTATION 


L 

s 

0.004650  H 

0.002199  H 

L 

0.000310  H 

m 

•■ff 

0.01553  H 

Sd 

0.007152  H 

‘■Dq 

■\f 

0.006293  H 

0.008334  H 

^ADD 

0.004716  H 

M.  * 

0.002678  H 

ADli 

^FDD 

0.008327  H 

*"al 

0.000367  H 

•-fl 

0.000774  H 

*"081 

0.000205  H 

*'Dql 

0.000278  H 

R 

3 

0.300  ohm 

0.2665  ohm 

'^Dd 

1 2 ohms 

R_ 

1 0 ohms 

Dq 

TABLE  9 

SATURATION  FACTOR  DATA  FOR  CONSTANTS  FROM  COMBINATION 
OF  TEST  AND  COMPUTATION 


where 

A = 0.012443137 
3 = 12.430776 

c = -38.595227 
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The  products  of  machine  constants  and  the  saturation  factors  for  test 
constants  and  for  the  Ai Research  computer  program  ROUND  constants^do  not 
greatly  differ.  The  transient  performance  derived  from  the  model  is  nearly 
the  same  whether  the  test  constants  or  the  AiResearch  aesign  program  ROUND 
constants  are  used. 
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SECTI ON  VIII 

DISCUSSION  OF  MODEL  SIMULATION 


Two  approaches  were  used  to  solve  the  set  of  differential  equations 
describing  the  six-circuit  model:  (1)  the  SCEPTRE  program  and  (2)  an 

Ai Research  program. 

Equation  (2)  can  be  written  as: 

..  dM, 


where 


and 


if.  _ df  d (MMF) 
dt  ~ d(MMF)  ’ dt 


d (MMF)  ^ ^^^d  ^ ^ ^ MMFq  _ d (MMFq) 

dt  MMF  ’ dt  MMF  ’ at 


d(MMF^)  di^  d^  ^ 

cit  dt  sfdt  2 


d i 


dt 


cos  Xt 


d i ^ ' c 

+ — ^ cos(jat  - 120)  + ^ cos  (-t  + 120) 


dt 

k 


af  r- 

- “ L' 


sin  x't  + i,  sin  (xt  - 120) 
a b 


+ i ^ sin  (xt  + 1 20) 


(47) 

(48) 


(43) 


d (MM Fq)  _ , 'dc  - *^af  ' ^ ' 
dt  ~ af  dt  2 
d i 

+ — - si  n (xt  - 1 20) 
dt 


i cos (xt+1 20) 


^ ' b 

sin  xt  + T7~  sin  (xt  - 1201 
at  ° 


^f 


I COS 

a 


xt  ' b cos  (xt  - 1 20) 
(50) 


Substituting  equations  (48)  through  (50)  into  equation  ^47) 


di^  di,  di  di,.  dip, 

= c.  .c^^  Co,—  -.o,l_Da  -C, 


— = c ~ ^ 

dt  ^a  dt 


dt 


(51 
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df 

d (MMF) 
MMF, 


MMF 


> ^3  F 

= f -j-  (P^  cos  out  - sin  u)t) 

' ^af  r 1 

C,  = f'  -:r-  P,  cos  (uat  - 120)  - P s in  (u)t  - 120) 

b 2 [ d q J 

^ r 1 

C = f'  — r-  P.  cos  (o)t  + 120)  - P sin  (tut  + 120) 


Sd ' f' 

Sc  = f' 


C = -f 

X 


I Pd  I i^  sin  uut  + i sin  (ujt  - 120) 

+ i sin  (ujt  T 120)  + P i cosxt+i , cos  (j 

c q [ 3 b 

+ cos  ('JUt  + 120) 


SCEPTRE  PROGRAM  APPROACH 


Since  — contains  derivatives  oF  state  variables  (i.e.,  — r-,  etc. 
^ ^ ^ ^ Q t 

term  — M i must  be  included  in  SCEPTRE  as  a Fictitious  vottaae: 
dt  s 


(L,  + fMs)  ^ + Rl  = (E 


- MS  ^ 1 - V) 


Iteration  and  numerical  integration  can  then  be  carried  out  like 
ordinary  SCEPTRE  problem. 

A I RESEARCH  PROGRAM  APPROACH 


I 


1 

L 
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where 


C Y, 
a 1 

'c'^1 

S'"! 

Sd^ 

Y, 
Dq  1 

C Y., 
a 2 

'b'^2 

'c'^2 

Sd'^2 

Sq^2 

C Y, 
a 3 

^c'^3 

Sd^3 

Sq^^3 

C Y, 
a 4 

C Y, 
c 4 

Sd\ 

Sq\ 

C Y, 
a 5 

^c^5 

'f^5 

''Dd^S 

Sq^S 

C Y, 
a 6 

'^b'^6 

'c'^6 

'^f^6 

Sd'^6 

Sq'^6  _ 

and  is  given  in  Equation  (60) 

Equat  ion  (46)  then  becomes 

(t,  - fM^  ^ C)  ^ = E - V - (R  + f ^ - C^\)  1 (65) 

A modified  Euler  technique  can  now  be  used  to  integrate  and  solve  this 
set  of  differential  equations.  The  matrix  C is  in  general  nonsymmetr ica 1 , 
and  is  added  to  the  symmetrical  inductance  matrix  (L^  + fMs) • The  resultant 
matrix  is  therefore  also  nonsymmet r i ca  I . This  is  the  basic  reason  the  SCEPTRE 
program  cannot  solve  Equation  (65). 

Since  no  circuit  change  is  allowed  in  SCEPTRE,  the  thyristor  switches  are 
simulated  by  resistors  that  change  in  value  depending  on  the  direction  of 
current  flow.  This  requires  extra  computer  time  to  iterate  into  current 
values.  It  also  results  in  leakage  currents  through  resistors  that  do  not 
actually  exist.  This  current  leakage  decreases  the  accuracy  of  the  solution. 

The  AiResearch  program  needs  no  iteration.  It  actually  simulates  the 
switching  phenomena  and  also  uses  much  less  computer  time  than  the  SCEPTRE 
program. 
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SECTION  IX 

COMPARISON  OF  TEST  AND  SIMULATION  RESULTS 

The  present  computer  program  'would  need  considerable  modification  in 
order  to  simulate  the  line-to-line  and  three-line  faults.  Therefore,  for 
load  cases  5 and  6,  double- 1 i ne-to-ground  fault  and  three-line-to-ground  fault 
cases  were  simulated.  The  tests  had  been  completed  on  the  experimental  alter- 
nator before  this  limitation  in  modeling  had  been  evaluated;  therefore,  j 
direct  test  comparison  for  the  model  results  was  made  in  load  case  5* 

Figures  29  through  35  are  computer  plots  of  transient  responses  for 
various  load  cases  using  the  AiResearch  mode  solution  program  ana  the  ROUND 
computer-program-calculated  generator  constants.  The  figures  are  discussed 
below. 

Figure  29 — Step  increase  of  applied  field  voltage.  The  corresponding 
test  result  is  shown  in  Figure  I7.  The  envelopes  of  the  transient 
excursions  in  Figure  I7  are  shown  dashed  in  Figure  29. 

In  the  following  figures,  dash  lines  shovy  the  envelopes  of  the 
test  transients  as  above. 

Figure  30 — Step  decrease  of  applied  field  voltage. 

Figure  31 --Normal  field  excitation,  sudcenly  applied  load.  The 
corresponding  test  result  is  shown  in  Figure  19.  In  the  test,  due 
to  misfiring  of  the  switch,  phase  C load  did  not  connect  until  two  | 

c'/cies  after  the  Phase  A oad  is  applied.  For  this  reason  t'e  tes: 
and  computer  results  do  not  exactly  match  during  a brief  transient 
period. 

Figure  32  — S ing le- 1 i ne- to-g round  fault. 

Figure  53--Doub le- 1 ine-to-ground  fault.  Test  results  are  not  shov/n 
for  comparison  in  this  case,  as  noted  above. 

Figure  34--T'nree- I i ne-to-qround  fault. 

Figure  3S--0ne  line  open. 

As  can  be  seen  from  the  figures,  the  transient  performance  precictec  by 
the  model  is  in  good  agreement  with  the  test  results. 
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Figure  36  shows  single-line-to-ground  fault  case  using  the  AiResearch 
model  program  and  the  combination  test  and  manually  computed  alternator  con- 
stants. As  this  comparison  shows,  these  results  are  very  close  to  those 
shown  in  Figure  32  (using  the  AiResearch  modal  program  and  the  AiResearch 
computer  ROUND  output  constants). 

The  SCEPTRE  program,  although  requiring  more  computer  time,  does  yield 

the  same  output  as  that  of  the  AiResearch  model  program.  Figures  37  through 

44  show  the  transient  responses  of  i^,  MMF , i,i,V  ,,V  .,i„,  and 

f ’ a c ca I ccl  Dd 

ipq,  respectively,  for  a single-line-to-ground  fault  using  the  AiResearch 
computer  ROUND  output  constants.  These  results  also  agree  closely  with  those 
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I iijuic  381).  Rcsultcint  MMr,  3 to  t>  Cycle  Run 
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Kiijure  ^Oc.  PhJse  t Cur  rent,  6 to  y Cycle  Kun 
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Kiijure  'tKi.  Phase  A Capacitor  1,  Voltaye,  0 to  3 Cycle  Run 


Filjuie  ^i2ci.  Plidse  C Copocitor  1 , Vollcige,  0 to  3 Cycle  Run 
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Figure  ^3c.  Direct  Damper  Current,  6 to  y Cycle  Run 
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SECTION  X 
CONCLUSIONS 

Results  of  the  alternator  test  and  computer  modeling  work  performed  were 
positive.  The  following  comments  and  conclusions  are  pertinent. 

The  test  alternator  and  load  i ng/contro  1 circuits  performed  as  intended 
with  no  serious  or  unforeseen  problems  in  construction  or  use.  Although 
precise  resonant  operation  was  not  realized,  the  alternator  was  operated  at 
400  Hz  with  tv/o  circuit  resonant  frequencies  so  that  half-cycle  current 
conduction  angles  of  approximately  125  and  160  electrical  degrees  were 
achieved.  (Note;  emphasis  in  the  program  was  on  the  test  and  modeling  cor- 
relation activities  in  order  to  verify  alternator  design  methods  and  approach. 

Achieving  precise  ac  resonant  operation  was  a secondary  object. ve.) 

Based  on  the  test  and  modeling  results,  the  AiResearch  alternator  design 
program  ROUND  does  develop  machine  saturation  data  and  constants  that  can  be 
utilized  to  confidently  predict  alternator  transient  behavior. 

Transient  simulation  results  were  checked  using  two  sets  of  alternator  j 

saturation  and  reactance  constant  data;  one  directly  taken  from  the  AiResearch  j 

ROUND  alternator  design  program,  and  one  based  on  actual  test  of  the  alternator  j 

supported  by  manual  calculation.  These  two  input  data  sets  generally  were  in 
close  quantitative  agreement  and  yielded  similar  results.  It  should  nonethe- 
less be  noted  that  the  accurate  quantitative  prediction  or  the  determination  . 

by  test  of  al^fmator  damper  circuit  constants  (particularly  for  nondiscrete 
elements  like  a solid  rotor  core)  are  presently  dependent  on  laraely  empirical  i 

i 

treatment.  This  is  an  area  of  possible  interest  for  further  analytical  inves- 
tigation. 1 

The  a I ternator/load  model  transient  behavior  cases  were  solved  using  two  | 

distinct  mathematical/computer  approaches;  the  SCEPTRE  program  and  an  | 

AiResearch  program.  Results  for  a given  input  data  set  were  in  very  close  j 

quantative  agreement.  .] 

While  SCEPTRE  is  capable  of  solving  the  al ternator/load  model  cases  and 
can  do  so  with  good  accuracy,  the  computing  time  is  still  excessive.  Accuracy 
of  SCEPTRE  results  is  essentially  tied  to  the  convergence  and  integration 
step  size  criteria  selected.  These,  in  turn,  directly  affect  computer  run 
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time  and  cost.  Even  when  rather  inexact  solution  is  allowed,  the  SCEPTRE 
program  requires  long  run  time.  Furthermore,  the  switching  phenomena  cannot 
be  exactly  treated.  Finally,  the  output  plot  data,  while  it  can  be  interpreted, 
is  not  convent ient ly  graphic  or  complete.  These  are  not  serious  defects,  but 
should  be  noted. 


